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Abstract. We study the interaction of low energy neutrinos on nuclei that 
spontaneously undergo beta decay showing that the product of the cross section times 
neutrino velocity takes values as high as 10 -42 cm 2 c for some specific nuclei that decay 
via allowed transitions. The absence of energy threshold and the value of the cross 
section single out these processes as a promising though very demanding approach for 
future experiments aimed at a direct detection of low energy neutrino backgrounds 
such as the cosmological relic neutrinos. 
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1. Introduction 

The interaction of an electron (anti)neutrino with a nucleus N naturally undergoing 
beta (positron) decay to the daughter nucleus N' 

( ^ + N-^N' + e ± , (1) 

shows the remarkable property of having no energy threshold on the value of the 
incoming neutrino energy. Indeed, in these cases the energy balance of the corresponding 
beta decay reactions is such that M(N) - M(N') = Qp > with M(N), M(N') the 
mass of neutral atoms. Neutrino interaction of this type is thus always energetically 
allowed, no matter the value of the incoming neutrino energy E v . Reaction (CQ) will be 
denoted in the following as Neutrino Capture on Beta decaying nuclei (NCB). In the 
limit of vanishing values of neutrino mass m v and E u the neutrino contributes to (JT]) 
uniquely via its lepton flavor quantum number and in this case the electron in the final 
state has exactly the beta decay endpoint energy Qp. However, for finite m u the electron 
kinetic energy is Qp + E v > Qp + m u , while electrons emerging from the analogous beta 
decay has at most an energy Qp — m u , neglecting nucleus recoil energy. A minimum 
gap of 2m u is thus present and this at least in principle allows to distinguish between 
beta decay and NCB interaction. 

In neutrino physics, NCB represents to date the only known reaction able to 
unambiguously detect electron (anti)neutrino having arbitrary low energies. Neutrinos 
with energies between 10 and 100 keV are still undetected even if some measurements 
have been proposed [2, El El E] • Low energy neutrino physics is indeed crucial to clarify 
our picture of fundamental particle interactions. Compelling results about neutrino 
mixing [5J and the corresponding non zero neutrino mass evidence raise a series of 
questions that could be addressed with the aid of the NCB process. The value of the 
absolute neutrino mass scale and the existence of low mass sterile neutrino mixing with 
ordinary flavors are examples of what can be done exploring this low energy regime. 

The idea of using NCB to measure the cosmological relic neutrino background 
predicted in the framework of the Hot Big Bang model was already advocated many 
years ago in [6J. The original idea was that if relic neutrinos have a large chemical 
potential /x, then both beta decays and NCB would show a depletion in the electron 
(positron) energy spectrum in a gap of order \x around the zero neutrino mass endpoint 
Qp. Presently, we know that these neutrinos have a number density of order n u ~ 50 
cm -3 neutrino (or antineutrino) per flavor and are characterized by a very small mean 
kinetic energy, of the order of 6.5 T^/m u or 3.15 T v for nonrelativistic and relativistic 
neutrinos, respectively with T„ = (4/H) 1 / 3 ^ ~ 1.7- 10~ 4 eV. Big Bang Nucleosynthesis 
constrains relic neutrino-antineutrino asymmetry, i.e. the ratio fi/T u for electron 
neutrino to be very small fi/T u < 0.1, see e.g. [TJ [8], and this bound also applies 
to all active neutrino flavors because of the oscillation mechanism [9]. This implies that, 
unless more exotic scenarios are considered, as a larger amount of relativistic degrees 
of freedom in the Early Universe, the effect of neutrino degeneracy in beta decays and 
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NCB is too small to be detected experimentally. However, as we mentioned already for 
massive neutrinos a similar gap around Qp is expected of the order of twice the neutrino 
mass, which for m v ~ 1 eV is several orders of magnitude larger than the corresponding 
effect due to neutrino-antineutrino asymmetry. 

Presently, neutrino mass in the eV range is still allowed by data. Indeed, oscillation 
experiments only provide a lower limit to the mass of (at least) one neutrino mass 
eigenstate of the order of 0.05 eV [10] . while direct measurements of electron energy 
spectrum in 3 H decay gives m u < 2 eV [H] H2] . A large improvement in this respect 
will be provided by the KATRIN experiment, whose expected sensitivity is 0.2 eV [13J. 
On the other hand, cosmological data from Cosmic Microwave Background anisotropics 
and Large Scale Structure power spectrum provide an independent bound on the sum 
of neutrino masses which depending on the particular model adopted and the number 
of free parameter lies in the range 0.3 -r- 2 eV, see e.g. [H] for a recent review. As we will 
argue in the following, if m v is in the eV range, future NCB experiments could represent 
an almost unique way to detect cosmological neutrinos. 

It is also worth mentioning that there are several possible sources of low energy 
neutrino fluxes which might be fruitfully studied using NCB. One example is provided 
by the thermal neutrino flux from the Sun, due to Compton pair production, free-bound 
electron transition and plasmon decay, which produce a neutrino/antineutrino flux in 
the energy range 1 eV - 10 KeV, with an integrated flux of the order of 10 6 — 10 7 cm' 2 s _1 
|15j . Neutrino diffused fluxes in the 10 KeV - lMeV energy range are instead expected 
from galactic and extra-galactic star nuclear burning [16J, as well as from POP III stars 
[PT] and for this latter case, their detection may represent a first direct probe of this yet 
unseen metal-free first generation of stars. 

Of course, the strength of neutrino interaction with ordinary matter is very small 
and NCB makes no exception. In order to design an experiment probing yet unexplored 
neutrino backgrounds a detailed estimate of the cross section is mandatory. We consider 
in the following a method to evaluate NCB cross section that makes use of beta 
decay observables in order to minimize the uncertainties due to nuclear matrix element 
evaluation. This procedure allows to obtain an accurate numerical estimate of the 
cross section in a large variety of cases. When possible, a comparison with previous 
calculations made for some specific nuclei by previous authors [HJ [19], [20] will also be 
performed, though NCB interaction in the low energy neutrino limit has not yet been 
systematically studied in literature. 

This paper is organized as follows. In Section 2 we recall beta decay formalism and 
derive the expression for the NCB cross section which is then evaluated for several cases 
in Section 3. In Section 4 the NCB to beta decay rate ratio is studied for the particular 
case of the cosmological relic neutrino background. Finally, in Section 5 we draw our 
concluding remarks. 
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2. Neutrino cross section on /3 ± decaying nuclei 

NCB and its corresponding beta decay are essentially the same phenomenon. In the case 
of unobserved polarization the two processes have the same invariant squared amplitude. 
In order to derive an expression for the NCB cross section we make use of the beta 
decay formalism, evaluating invariant amplitudes using the description of the fermions 
in terms of spherical waves and the computation of the corresponding nuclear state 
transition matrix elements. This approach is based on the so called long wavelength 
limit approximation p v R <C 1, where p v is the momentum of the incoming neutrino and 
R the nuclear radius. Due to the typical values of R, this approximation holds for NCB 
if the energy of the incoming neutrino satisfies the relation E v < 10 MeV. 

Assuming an isotropic neutrino flux, the NCB integrated rate is given by (hereafter 
we use natural units) 

/(Ppi/ 
o-kcbVv f(p u ) , (2) 

where cr NCB is the NCB cross section, v u and p v the modulus of neutrino velocity and 
momentum, respectively and f(p v ) the distribution function. For example, in the 
particular case of cosmological relic neutrinos it keeps a standard Fermi-Dirac form, 
f(p u ) = [exp(p v /T u ) + up to very small non thermal corrections [21] . Following the 
formalism of [22], for spin averaged initial state and unobserved polarization we have 

G 2 

o-kcbVv = —p e E e F(Z, E e )C(E e , p v ) v , (3) 

7T 

where F(Z, E e ) is the Fermi function. With E e and p e we denote the energy and 
momentum of the outgoing electron respectively, where 

E e = E u + Qp + m e = E v + m v + W Q , (4) 

with W the corresponding beta decay endpoint. Finally, the rate can be expressed as 
an integral over the electron energy 

Q2 roo 



K = tt^ r PeE e F(Z, E e )C(E e , p v 

^ JWo+1m v 



'W +2m v 

■ E u Pu f \Pu) dE e , (5) 

where Eq. (jl]) should be used to express neutrino energy and momentum as function of 
E e . 

The NCB rate contains the nuclear shape factor C(E e ,p 1/ ) u , an angular momentum 
weighted average of nuclear state transition amplitudes, which depends upon the nuclear 
properties of the parent and daughter nuclei, thus involving the calculation of nuclear 
matrix elements. Its general expression is, see [22] for further details 



C(E e ,p u )p= \ ke Mx(k e ,k u ) + m 2 K (k e ,k UJ 

f^e feu •jK 



k e E e 



M 2 K (k e ,k u )m 2 K (k e ,k u ) , (6) 
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where k e and k v are the electron and neutrino radial wave function indexes (k = j + 1/2), 
K represents the nuclear transition multipolarity (\k e — k u \ ^ K ^ k e + k v ) and 
M^ike, k u ), m 2 K {k ei k u ) are nuclear form factor functions, which are typically evaluated 
using N-body nuclear transition matrix elements in the impulse approximation. This 
represents the main source of uncertainty for <7 NCB . 

On the other hand, NCB rate is strongly related to the corresponding beta decay 
process, whose rate is given by 
G 2 r w ° 

X P = t4 / PeE e F(Z } E e )C(E e} p v )pE u p v dE e} (7) 

a simple relation holding between the beta decay and the NCB shape factors [22] 

C(E e ,p v ) v = C(E e ,-p v ) fi , (8) 

though both variables have different kinematical domains in the two processes. 

The beta decay rate provides a relation expressing the mean shape factor, defined 

as 

_ 1 r w ° 

C p = - / p e E e F(Z, E e )C(E e ,p v ) p E v p v dE e , (9) 

J J m e 

in terms of observable quantities, W Q and the half-life t%/2 

27r 3 ln2 . , 

/*i/a = 7^-, (io) 



where / is defined as 

-Wo 



rw 

f= F(Z,E e )p e E e E uPl/ dE e . (11) 

Jm e 

We therefore obtain 

<W„ = 27r 2 ln2 p e E e F(Z,E e ) C ( Ee ' P " )u , (12) 

where the two factors C(E e , p u ) u and C p depend upon the same nuclear transition matrix 
elements. 

It is useful to define the quantity A as 

a = L£i HQ) 

p e E e F(Z,E e )C(E e ,p u )^ 1 ; 

and note that it contains the ratio of NCB and beta decay shape factors 

A _ [ W ° C(E' e ,p' v )p p> e E' e F(E' e , Z) 

J me C(E e , Pu ) u p e E e F(E e ,Z) W ^> li4j 

where a prime denotes all variables depending on E' e which should be integrated over. 
The A is a function of E v only, once the target nucleus characterized by Qp and Z is 
given. The NCB cross section times neutrino velocity can then be conveniently written 

as 

27r 2 hi2 

VkcbVv = -7—. , (15) 

A ■ ti/2 
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with the product ftyz characterizing the various beta decay processes. Values of 
log(/ti/ 2 ) are reported e.g. in [23J. 

As we will see, in some relevant cases the evaluation of A is particularly simple so 
that Eq. (fl5|) can be computed in an exact way. In all cases where this is not possible, 
systematic uncertainties affecting the nuclear matrix element evaluation are expected 
to largely cancel in the shape factor ratio appearing in A, and this suggests that the 
use of expression ffl5l) might be very useful to get a reliable estimate of the NCB cross 
section. 



2.1. Superallowed transitions 

Beta decays in which the nucleus undergo a transition taking place between members 
of the same isospin multiplet are called superallowed. Due to the large superposition 
between initial and final nuclear states these decays have the lowest known /ti/2 values. 
In particular, if the transition is of the + — > + type only the vector current contributes 
to the decay. In the limit of CVC hypothesis the vector form factor is given by the Fermi 
matrix element and the shape factor can be written as 

C(E e , Pv ) = = (F) 2 = (T — T 3 )(T + T 3 + 1) , (16) 

where T and T 3 are the isospin quantum numbers. 

In all other transitions happening inside the same isospin multiplet (J 71 — > J 71 " , J ^ 
0) both the vector and axial form factors contribute and in this case the shape factor 
can be written as 

C{E eiPv ) = TO 2 + |^)| 2 = (F) 2 + (W^) 2 (GT> 2 , (17) 

where g Ay is the axial (vector) coupling constant. As first approximation both form 
factors can be evaluated using the initial and final state spin and isospin wavefunctions. 
However, a precise calculation requires the evaluation of higher order correction to 
nuclear matrix elements. 

We notice that for all superallowed decay the form factors do not depend on E e . 
In this case A can be written as 

A = PeEj(Z, E e ) ' (18) 

and we get for NCB processes 

9l p e E e F(Z,E e ) , , 

a nCB v u = 2vr 2 In 1 Pe e f ) e) , (19) 

7*1/2 

which only depends on the half-life and the Qp of the corresponding beta decay. 

To illustrate the difference between expressions Q and f|T5|) we consider the specific 
case of Tritium which decays to 3 He with Q@ = 18.591(1) keV [24J and an half-life of 
t 1/2 = 12.32(4) years [25]. We adopt (F) 2 = 0.9987 and (GT) = ^3 • (0.964 ± 0.016) 
from [26] assuming a total 1.6% systematic uncertainty on the Gamow- Teller matrix 
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element evaluation. Using expression (j3j) and ffTTjl with Gp = 1.16637(1) x 10 5 GeV 2 , 
g A = 1.2695(29) and \V ud \ = 0.97377(27) Q [5J we find 

a NCB ( 3 H) — = ( 7.7 ± 0.2 ) x 10" 45 cm 2 , (20) 

c 

in case of massless neutrinos and in the limit of p v — > 0. We notice here that the 
main contribution to the quoted error is due to the uncertainty on the evaluation 
of the Gamow- Teller matrix element. This result is in agreement with previous 
calculations [IHlETj. On the other hand, evaluating <7 NCB ( 3 H) using expression ( [191) 
gives 

a NCB ( 3 H)— = ( 7.84 ± 0.03 ) x 10~ 45 cm 2 , (21) 
c 

where the error is entirely due to the experimental uncertainties on Qp and ti/ 2 and the 
effect of matrix element uncertainty is now fully removed. 



2.2. Allowed transitions 

The neutrino cross section for nuclei that undergo allowed beta decay can be evaluated 
only using a first order approximation since the beta decay shape factor is written in 
general as 

C(E e ,p u ) = TO 2 + \ A F^\ 2 + 0(p e R)0(aZ) . (22) 

If only the leading terms are taken into account one gets the following approximation 

C(E e ,p u ) l3 = C(E e ,p u ) u = constant , (23) 

and from ([3]) one derives an order of magnitude estimate for the NCB interaction 
cross section. The effect of the higher order terms is however, not negligible and more 
precise estimates require nuclear matrix element evaluation for each specific transition. 
Nevertheless, we expect that at a higher level of accuracy with respect to the leading 
approximation as above can be obtained using ffl5|) and the fact that for the shape factor 
ratio is reasonable to assume 



C(E e ,p v ) p _ 1 
C{E e ,p v ) v 

Within this approximation A is obtained using expression ([TBI . 



(24) 



2.3. Unique K-th forbidden transitions 

Neutrino cross section on nuclei that decay via K-unique forbidden beta decay takes a 
simple form since there is only one nuclear form factor involved in the shape factor. In 
particular, its expression can be written as 

C(E e ,p v ) p =\ A F%_ u \ 2 

L 

X YsSl^PeRY^tPuR) 2 ^ , (25) 
n=l 

| We recall here that errors on g A and V u d are strongly correlated. 
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where K is the degree of forbiddeness of the decay, L = K + 1 and £>£ and A n are 
respectively, numerical coefficients and functions that can be evaluated as in [22J. In 
particular, we can define the following functions Ui for i-forbidden decays 

Ui(p e ,Pu) =Pl + A 2 Pe> ( 26 ) 

U 2 (Pe, Pu) =Pt + y ^Plpl + A 3 Pe , ( 27 ) 

u 3 (Pe, P.) =pl + 7\ 2 ptp 2 e + 7\ 3 plpt + ^Pl ■ (28) 
The corresponding shape factors are then given by 

Rj 



C(E e ,p u ) 



_(2i + 1) 



•af(q) 
I c»+i)»i| 



Ui(Pe,Pv 



(29) 



Since there are only even powers of the neutrino momentum p v these expressions are 
valid both for beta decay and NCB interaction. The quantity A can be written using 
these expressions as 

^u t (p' el p'ME' e F(Z,E' e ) 



Ai 



Kp'AK 



u i (p e ,p v )p e E e F(Z,E l 
and results to be independent of nuclear form factors 



(30) 



3. Estimating the cross section 

We have evaluated NCB cross section for several nuclei using parametrization of the 
Fermi function and of the radial wave function coefficients X n in the shape factor as in 
[22J and [28] to account for finite nuclear size effect. Two numerical programs have been 
used to evaluate numerically F(Z,E e ) and A n , implemented with Mathematica [29] 
and C++ in a Linux environment, using a high precision calculation of the Gamma 
function [30]. In both cases the algorithms have been checked against each other and 
compared with the values reported in [31] in the range where these are available. A 
total accuracy of 10~ 3 was achieved on F(Z, E e ) and A n for values of Z up to Z ~ 40, 
while for higher values of Z and at large values of electron (positron) energy (E e > 10 
MeV) the accuracy worsen to a few percent. 

To illustrate the behavior of a NCB v u as a function of the nucleus we first consider 
in Figure [1] the ratio Q\jA for superallowed and unique-forbidden transitions in the 
case of massless neutrino and in the limit of small neutrino momentum. Notice that as 
expected, this ratio is only weakly depending on Q@ as from ([141) A grows as Q%. 

As we mentioned in Section II, the /ti/2 value characterizes the particular beta 
decay process corresponding to the NCB. Values of log(/ii/2) are reported in [23] with 
typical mean values of 3.44 for superallowed transitions, 5.5 for allowed transitions, 9.5 
for first unique forbidden transitions, 15.6 for second unique forbidden transitions and 
finally, 21.1 for third unique forbidden transitions. Using these results the NCB cross 
section can be derived using (115j) once the value of / is known. We show the typical 
values of this parameter in Figure [2] for several choices of Z. 
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Figure 1. Values of 1/A ■ versus Qp evaluated for (3 decaying nuclei in the 
limit of p v — ► 0. The four curves represent from bottom to top superallowed, first 
unique forbidden, second unique forbidden and third unique forbidden transitions, 
respectively. Curves are shown for Z — 20 and a nuclear radius given by R — 1.2A 1 / 3 
fm, where A — 2.5Z 



The behavior of a" NCB tv versus the incoming neutrino energy is shown in Figures [3] 
and H] for (3~ and (3 + decaying nuclei, respectively and for a particular choice of Z. 
Notice that it reaches a plateau at low neutrino energies, whose value strongly depends 
upon the kind of transition (super- allowed, allowed, etc.) and Qp. We have performed 
an extensive calculation of the NCB cross section for all beta decaying transitions listed 
in the ENSDF database [32]. A total of 14543 decays have been analyzed, 6288 (3~ and 
8255 electron capture and (3 + . We restricted our attention to both allowed and unique 
forbidden decays having branching ratios greater than 5%, namely 1272 f3~ decays 
and 799 f3 + decays. The corresponding a KCB v v for low incoming neutrino momentum 
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Figure 2. Values of / • 1/Q| for (3 decaying nuclei in the limit of p u — ► 0. The curves 
are evaluated using different values of the parent nucleus Z, from bottom to top (top 
to bottom in case of /3 + ) Z=20,40,60,80, respectively. The nuclear radius is taken to 
be R = 1.2A 1 / 3 fm, where A = 2.5Z 



are shown in Figure [5j There are several nuclei spanning a wide range in for 
which interesting high values are reached in the range 10~ 41 — 10~ 43 cm 2 . Some of 
these are reported in Table [U were we collect the results for the best know + — ► + 
superallowed decays. As explained in section II, our estimate is based on the knowledge 
of experimental values of Qp and £1/2 which in this case are both known with a great 
accuracy as this class of decays is used to test the CVC hypothesis within the Standard 
Model [33]. 

As we will discuss in details in the next Section, any use of NCB processes to study 
low energy neutrino fluxes is crucially related to the issue of rejection of background 
events represented by the corresponding beta decay process. For a given incident 
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Figure 3. The product Uncb^i; f° r &~ decaying nuclei versus neutrino energy. 
Typical values for log(/£) values have been assumed [23] as follows: (a) allowed, 
log(ft 1 / 2 ) = 5.5, (b) first unique forbidden, log(/i 1 / 2 ) = 9.5, (c) second unique 
forbidden, log(/ti/ 2 ) = 15.6, (d) third unique forbidden, log(/i 1 / 2 ) = 21.1. The 
three curves refer to different Q^-values, solid line for Qp = 10 -3 MeV, dashed line 
for Qj3 = 1CP 1 MeV, dotted line for Qp = 10 MeV. Curves are for Z = 21 and nuclear 
radius given by R = 1.2A 1 / 3 fm, where A = 2.5Z 

neutrino flux, the ratio of the NCB to decay events is proportional to a^ c ^{y v /c) ■ t±/2 , 
so that nuclei with the highest value for this combination might give the best chances 
if used in a future low energy neutrino detection experiment. We report these nuclei in 
Table El Among these, we recognize two isotopes, namely 3 H and 187 Re, widely used in 
the past and ongoing experiments or representing a future possibility for calorimetric 
neutrino mass experiment [3l] . 
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Figure 4. The behavior of <j kc -bVu for (3 + decaying nuclei versus neutrino energy. We 
assume the following values of log(/ii/ 2 ) [23]: (a) superallowed, log(/ii/ 2 ) = 3.44, (b) 
first unique forbidden, log(/i 1 / 2 ) = 9.5 (c) second unique forbidden, log(/ti/ 2 ) = 15.6, 
(d) third unique forbidden, log(/ii/ 2 ) = 21.1. The three curves refer to different Qp- 
values, solid line for Q — 10 -3 MeV, dashed line for Qp — 10 _1 MeV, dotted line for 
Qp = 10 MeV. Curves are for Z = 21 and nuclear radius given by R = 1.2A 1 / 3 fm, 
where A = 2.5Z 



4. NCB versus (3 decay: the case of cosmological relic neutrinos 

We now consider a possible application of NCB to detection of the background of 
cosmological relic neutrinos. Actually, this represents one of the most ambitious 
challenges in modern cosmology. There are two intertwined issues which should be 
discussed, i.e. the event to background rate and the energy resolution. First of all, 
though as we stressed several times the NCB process is with no energy threshold, 
nevertheless, the ratio of NCB event number to corresponding beta decay events is 
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Table 1. The product cr NCB (iv/c) for the best known superallowed + — * + 
transitions. Numerical values for Qp and partial half-lifcs are taken from 33J . The 
value of / is calculated adopting the parametrization of the Fermi function of [28] . 



Isotope 


/O l~\ r ^\T\ 

Qp (keVJ 


riali-Iiie (sec) 


Cncb(^/c) ' 


io c 


885.87 


1320.99 


5.36 x 10~ 3 


14Q 


1891.8 


71.152 


1.49 x 10~ 2 


26m^y 


3210.55 


6.3502 


3.54 x lO- 2 


34 C1 


4469.78 


1.5280 


5.90 x 10" 2 


38m j£ 


5022.4 


0.92512 


7.03 x 10~ 2 


42 Sc 


5403.63 


0.68143 


7.76 x 10~ 2 


46y 


6028.71 


0.42299 


9.17 x 10~ 2 


50 Mn 


6610.43 


0.28371 


1.05 x lO" 1 


54 Co 


7220.6 


0.19350 


1.20 x 10- 1 



Table 2. Beta decaying nuclei that present the largest product of jncb^/c) ■ t 1 / 2 for 
low neutrino momentum and have a /3 decay branching fraction larger than 80%. 



Isotope 


Decay 


Qp (keV) 


Half-life (sec) 


<Txob(v v /c) (10 41 cm 2 ) 


3 H 


P~ 


18.591 


3.8878 x 10 8 


7.84 x 10" 4 


63 Ni 


P~ 


66.945 


3.1588 x 10 9 


1.38 x 10~ 6 


93 Zr 


Pr 


60.63 


4.952 x 10 13 


2.39 x 10~ 10 


106 Ru 


Pr 


39.4 


3.2278 x 10 7 


5.88 x 10" 4 


107p d 


P~ 


33 


2.0512 x 10 14 


2.58 x 10~ 10 


187 Re 


Pr 


2.64 


1.3727 x 10 18 


4.32 x lO" 11 


n c 


P + 


960.2 


1.226 x 10 3 


4.66 x 10~ 3 


13 N 


P + 


1198.5 


5.99 x 10 2 


5.3 x 10~ 3 


15Q 


P + 


1732 


1.224 x 10 2 


9.75 x 10~ 3 


18p 


P + 


633.5 


6.809 x 10 3 


2.63 x 10~ 3 


22 Na 


P + 


545.6 


9.07 x 10 7 


3.04 x 10~ 7 


45 Ti 


P + 


1040.4 


1.307 x 10 4 


3.87 x 10" 4 



typically very small. Using (TSJ) and ( 1151) we find 

M' im ^.4^(l»n^W (31) 
\p \p,-*o J In 2 \p,-*o A J 

where we have used the fact that relic neutrinos have a very small mean momentum of 

order T u with a spread of the same order of magnitude, and the fact that the product of 

NCB cross section times neutrino velocity gets an asymptotic constant value for small 

neutrino energies, see Figures [3] and HI In the case of 3 H we get 

A^( 3 H) = 0.66 ■ 10~ 2 %( 3 H) . (32) 
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Figure 5. <7 NCB (viy/c) for (3^ nuclei as a function of the decay Q^-value and evaluated 
in the limit of low p v . A total of 1272 f3~ and 799 (3 + decays having strong (3 branching 
ratio (BR > 5%) is reported for both allowed and forbidden decays 



A numerical estimate of the relative event rate of low energy neutrino NCB interaction 
with respect to beta decays can be obtained by using the results of Figured] in case of 
allowed and unique forbidden transitions. 

Despite of this disappointing result, at least in principle the experimental signature 
of NCB events is unambiguous as the electron (positron) in the final state has a kinetic 
energy at least 2m u above the beta decay endpoint energy. However, the finite energy 
resolution of any experimental apparatus and the extremely low cross section make 
relic neutrino detection via NCB a real challenge due to the large background events. 
In particular, for low neutrino masses smaller than the typical experimental energy 
resolution, it is really impossible to disentangle the few expected NCB events from the 
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large background of standard beta events. In this case NCB processes are of no use. 

On the other hand, in more optimistic scenario with comparable values of neutrino 
masses and experimental energy resolution, the situation could be much more promising. 
As an example, we consider a future experiment reaching an energy resolution A, and 
neutrino masses in the eV range. From (j2J) and (J7J), the ratio of the event rate A y g(A) 
for the last beta decay electron energy bin W a — A < E e < W D , compared with the total 
NCB event rate can be easily calculated, giving 

X^A) =9 2 C{3) ^k) (1 + 2^/A) 3 / 2 ' (33) 

where we have used that n u = 3£(3)T^/(47r 2 ) and the fact that Qp » A. We 
have checked that this expression is accurate at percent level for nuclei that undergo 
both allowed and unique-forbidden transitions with endpoint energy in the range 
1(T 3 < Q/3 < 10 MeV. This gives for example, the value A„/A/j(A) ~ 2.2 • 1(T 10 for 
A = 0.2 eV and m v = 0.5 eV. 

To get an estimate of the signal to background ratio we can proceed as follows. For 
A < m v , the expected background electron events which are produced by beta decay, 
yet having an energy which corresponds to the relic neutrino capture energy bin centered 
at E e = W + 2m u is suppressed by the exponential factor 

1 /■2m„/A+l/2 

p = -== / e~ x2/2 dx . (34) 

V27T j2mu/A-l/2 



Thus, imposing a signal to background ratio larger than unity corresponds to the 
condition 

HS V + wAr, - 1 ' (35) 

and this for a given neutrino mass provides the energy resolution which is necessary to 
achieve. A signal to noise ratio of order 3 is for example obtained if A = 0.2 eV for 
m v = 0.7 eV, while a smaller neutrino mass of 0.3 eV requires A = 0.1 eV. In these 
cases a total event number of order 10 is needed to get a 5-cr discovery claim. Presently, 
this energy resolution seems very hard to get. Nevertheless, if a large neutrino mass will 
be found by ongoing beta decay experiments such as KATRIN, it is not inconceivable 
that a future generation of experiments might reach energy resolution as low as 0.1 eV. 

Finally, we estimate the order of magnitude of the mass of detector required to 
see neutrino events from the cosmo logical background using NCB. For a mass M[g] 
expressed in grams, the expected total event rate is 

A„^M, (36) 

where Na is the Avogadro number and A is the atomic number of the decaying nucleus. 
Inserting numerical values we get for the molar rate 

2 g5 . 1Q -2 wWc , ^ 
10 4b cm 2 
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Table 3. The number of NCB events per year for 100 g of 3 H, taking into account 
the effect of gravitational clustering in the neighborhood of the earth, compared to the 
case of a standard homogenous Fcrmi-Dirac distribution with T v = 1.1 ■ 10 -4 eV (FD). 
We show for some value of neutrino mass the results for a Navarro, Frenk and White 
profile (NFW) and for present day mass distribution of the Milky Way (MW), using 
the local neutrino densities computed in [36] . 



m v (eV) 


FD (events yr x ) 


NFW (events yr" 1 ) 


MW (events yrs 1 ) 


0.6 


7.5 


90 


150 


0.3 


7.5 


23 


33 


0.15 


7.5 


10 


12 



As an interesting example, we consider the case of 3 H. From (1371) and using the results 
of Table [2J we estimate 7.5 events per year of data taking for a mass of 100 g. On the 
other hand, a very small result is obtained in the case of 187 Re, due to the tiny NCB 
cross section, o" NCB ( 187 Re)f iy /c ~ 10 -52 cm 2 . 

We conclude this Section by observing that all our findings are obtained assuming 
a standard and homogeneous relic neutrino background. However, massive neutrino 
density could be locally larger because of gravitational clustering. This effect in a Cold 
Dark Matter Halo could be relevant for order eV neutrino masses [35l [36]. Indeed, 
the neutrino density in the neighborhood of the earth is enhanced by a factor as large 
as 10-i-20 for m v = 0.6 eV or 3-^4 if m u = 0.3 eV [36] • On the other hand, neutrino 
distribution in phase space is very close to the standard homogeneous Fermi-Dirac result 
for neutrino masses smaller than 0.1 eV. To quantify the effect of gravitational clustering 
on NCB rate we again consider the expected event rate for 100 g mass of 3 H. Results 
are shown in Table [3j The enhancement is simply due to the larger integrated neutrino 
density n u . Actually, even for large neutrino masses the neutrino distribution has a 
mean momentum of the order of 5-=-8 10~ 4 eV and decreases exponentially for momenta 
larger than 10 -3 eV. This holds for both a Navarro, Frenk and White and present 
day Milky Way mass profiles, see Figure 6 in [36]. Therefore, the event rate can be 
always computed as in (I3TT) in the limit of vanishing neutrino momentum and taking 
into account the local higher neutrino density n u . This results in a remarkable increase 
of the expected event number for large neutrino mass and thus, for a fixed ratio A/m u , 
in a higher statistical significance due to the larger signal to (beta decay) background 
ratio. 

5. Conclusions 

The detection of low energy neutrino backgrounds, primarily those emitted at the freeze 
out of weak interactions in the early stages of the evolution of the Universe, is still beyond 
our present experimental capability. Nevertheless, there is a continuous struggle both 
in the theoretical and experimental physics communities in searching for new ways and 
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methods to detect these fleeting fluxes. 

In this paper we have reported a careful analysis of neutrino capture on beta 
decaying nuclei. These processes have the remarkable property of having no energy 
threshold on the incoming neutrino energy and thus they might represent a good and 
numerous class of interactions suitable for low energy neutrino detection. This idea was 
already suggested several years ago by S. Weinberg in [6], though in the framework of a 
degenerate massless relic neutrino background. Our study has been inspired by this and 
prompted by the fact that we presently know that neutrinos are massive particles and 
might have masses m v as large as a fraction of eV. Indeed, there is a gap of 2m u in the 
electron or positron energy spectrum separating the (few) events induced by e.g. relic 
neutrinos, from the large background of standard beta decays. At least in principle, this 
allows to disentangle the two processes. Of course, any possibility to translate this idea 
into a real experimental technique depends upon two crucial issues, namely the expected 
order of magnitude of NCB event rate as well as the required energy resolution of the 
outgoing electron (positron). 

The first issue has been considered in details in this paper. We have studied the low 
momentum limit of the NCB cross sections for more than two thousands beta decaying 
nuclei with a branching factor larger than 5%. Whenever possible, we have tried to 
reduce the effect of the uncertainty on the nuclear matrix elements by linking the NCB 
rate to the experimentally known properties of the corresponding beta decay, the half- 
life and the Qp value. Depending on the involved nuclear transition, the cross section 
times neutrino velocity in the low neutrino momentum regime spans several order of 
magnitude and interestingly, can be as large as 10~ 42 4- 10~ 43 cm 2 c for super-allowed 
decays. This means that the event rate can be remarkably large. For example, we 
found that with 100 g of 3 H one expects order 10 events per year due to scattering 
of cosmological relic neutrinos. Similar rates are also found for many super-allowed 
transitions. Interestingly, the event rate is expected to be even larger for neutrino 
masses of 0.3 4 0.7 eV due to the effect of gravitational clustering in the neighborhood 
of the earth. In this case the expected signal rate is in the range 204150 per year per 
100 g of 3 H. 

The second problem is the possibility to have a reasonable rejection of the 
background due to standard beta decay electrons/positrons. This depends on the 
incoming neutrino energy range which is under study and the experimental energy 
resolution. In the most demanding case of the extremely low momentum cosmological 
relic neutrinos, it is necessary to reach a sensitivity which is better than the value of 
neutrino mass. This can be foreseen as a plausible perspective for future experiments 
only if neutrinos have masses of order eV, thus in the so called degenerate scheme for 
neutrino masses, which is still allowed by all present data, though slightly disfavored by 
cosmological observations. For example, for m v ~ 0.5 eV and an energy resolution of 
0.1 4 0.2 eV a reasonable event to background discrimination is possible. For smaller 
neutrino masses it seems very hard if not impossible at all to use NCB processes for 
relic neutrinos. 
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Of course, further improvements of the background rejection could be achieved by 
a careful and more complete reconstruction of the kinematics of the observed events, 
including possibly the daughter nucleus recoil energy. As a further example, for polarized 
nuclei one could also measure the beta-gamma angular correlation in case the daughter 
nuclei de-excite via gamma emission. However, these detailed issues are beyond the 
aim of the present paper and might deserve further insights depending on forthcoming 
information about the absolute neutrino mass scale. 
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